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Abstract. We show that the tritronquee solution of the Painleve equation Pt, y" = 6j/ 2 + 2 
which is analytic for large z with arg 2 £ (— t^t) is pole- free in a region containing the 
full sector {2 7^ 0,arg2 6 [— ^?,7r]} and the disk {2 : \z\ < |^}- This proves in particular the 
Dubrovin conjecture, an open problem in the theory of Painleve transcendents. The method, 
building on a technique developed in [4], is general and constructive. As a byproduct, we obtain 
the value of the tritronquee and its derivative at zero within less than 1/100 rigorous error 
■ bounds. 

o 

_ 1. Introduction and Main result 

Understanding the global behavior in C of the tritronquee solutions (see below) of the Painleve 
equation Pi is essential in a number of problems such as the critical behavior in the NLS/Toda 
. lattices ([BJ, [7]) and the analysis of the cubic oscillator ([IS])- Considerations related to the be- 

havior of NLS/Toda solutions corroborated by numerical evidence led to the Dubrovin conjecture, 
see [6J (cf. also [2] and [7]). This conjecture states that the tritronquee solutions are analytic 
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in a neighborhood of the origin O and in a sector of width 87r/5 containing O. A number of 
partial results on this question have been obtained so far, see e.g. [8], [E]-[l4] but, in spite of the 
existence of an underlying Riemann-Hilbert representation, at the time of the present paper the 
conjecture is still open. 

The purpose of the present paper is to prove the Dubrovin conjecture alongside other results 
about the tritronquees. 

^ , We write Pi in the form 

0\ '. (1) y" = 6y 2 + z 

■ Tritronquees. There are exactly five solutions of ([T]) which are analytic for large z in a sector of 
width 87r/5. These special solutions called tritronquees are obtained from each-other through the 
five-fold symmetry of Pi, y i-> e ^/5 y (^ife/s^ } k = 0,...,4 (cf. [9], [11], M)- To understand 
their properties it is clearly enough to analyze one of them; we choose the solution y t uniquely 

■ defined by the property y(z) ~ i^z/6[l + 0(z~ 5 / 2 )] as z — > +00 (cf. e.g. Proposition 2 and 
Theorem 3 in [8]). The main result of this paper is 

Theorem 1. The tritronquee yt is analytic in the region 
(2) {z^O: argze [-3tt/5, tt] } U {z : \z\ < § } 

Corollary 1. Dubrovin's conjecture, stating that y t is free of singularities in {z : argz £ [— 1 7r, 7r] }, 
holds. 

Note 1. (i) The proof of Theorem [T] is organized as follows. We first establish the analyticity of 
y t in the region {z : \z\ > 1.7, arg z <G (— 37r/5, 7r)}, and continuity in its closure. We determine 
y t (zo) , y' t (zo) , zq = 1.7e 47r / 5 with errors smaller than 6.5 • 10~ 3 . We then show that any solution 
of Pi with y t (zo), yj(zo) hi this range matches a Maclaurin series solution of ([T]) with a radius of 
convergence bounded below by 37/20 = 1.85. The result follows. 

(ii) As a byproduct we find within rigorous errors bounds smaller than 1/100 the values of 
?/(0) and j/'(0), which arc also needed in applications, and these values are consistent with the 
numerical calculation by Joshi and Kitaev [8]. 

(ii) The lower bound 1.85 is not optimal, yet not far from the numerically obtained radius of 
analyticity, ~ 2.38 ([8]). The methods we use can, in principle, be adapted to rigorously calculate 
the tritronquee and the position of its first pole with any prescribed accuracy. 

We will also use a standard normalization of Pi [TJ, similar to the Boutroux form. After the 
change of variables 

' n5/4 / n . I Z I _ 



(3) * = ~W^' = Ve ( 1 -2^ + h ^ 
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Pi becomes 

ui l u 7 h 2 392 
(4) h" + -h! -h = — 



x 2 625cc 4 

Using the change of variables leading to (j4|), the aforementioned reference [8] implies that there is 
a unique solution ft. of (J4j) with the behavior h(x) = O (x~ 4 ) as x — > +ioo. 

2. The tritronquee y t on the anti-Stokes line z = e 47r / 5 M + 

First, we consider behavior of y t on the anti-Stokes line z = e 47T / 5 M + for \z\ ^ (30) 4 ^ 5 (= 
0.633 • • • ). We find yt(zo) and y' t (zo) for z = 1.7e i7r / 5 , accurately and with rigorous error bounds, 
later used to show that any solution of Pi within the proved range has y t with a power series at 
the origin with radius of convergence exceeding 1.85. After the substitution 

(5) h(x) =: — =-, 

y x 

d4| implies 

„ / 1 \ H 2 392 

(6) ff"- 1--, 5 = 



2 



4x 2 / 2^i 625x 7 / 
We define 

(7) 7 := {.-e G iIR + , |a;| > p} ; where p > 1 
and the Banach space 

(8) 5 7 = {J?GC(n j! C):||F||<oo} 
where || • || is the weighted norm 



(9) \\H\\ = sup 

xEQi 

We invert ([6]) in 5/ to obtain the integral equation 



, b ' 2 H 



(10) H{x) = H Q (x) + I wak(x - t) <j + ^-M } dt =: N [H] (x 



H{x) = H Q {x) + [ sinh(a; - t) [- 

</ zoo L 



4i 2 2^i 



where 

392 

(11) H (x)= sinh^-t) __d* 

^/ 200 

In the process of inverting no nonzero linear combination C\e x + C<x€T x may be added to the 

24 



625i 7 / 2 
lation C] 

right side of ([TOf as e ±:r ^ 5/. We note that |x| = p = 1 corresponds to |z| = 30 2 . / = 0.6331 • • 



Lemma 2. There exists a unique solution H to UU\) and therefore of (0|) in Sj. 

Proof. Using |sinh(a; — t)\ ^ 1 for x,t € Hi and (fTTj) we obtain 

784 784 
(12) \H (x)\^ — ^ ^||g |K 



3125|a;| 5 / 2 " " 3125 
Since \H(t)\ < |i|~ 5 / 2 it follows that 

\H\\ . Ilffl 



W[H]\\ < ll#o 



14|x| 7 / 2 9|x| 9 / 2: 



Therefore, it follows that 

llffll ll^ll 2 



14p 9P 2 ' 
and similarly, 

W Wi] - M [H 2 ] ||< + l|gl|l + 2 l|ff211 ) ||^ - ff 2 
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Consider now the ball B C Sj of radius (1 + e)||ifo|| for some e > to be chosen shortly. For any 
H , Hi , H2 G B we obtain from the inequalities above 



(13) 



||^[tf]K||#o|| 1 



i+f (1 

14p 



v 



{ii + 9^ 1+ ^' 



(14) \\M[H X ]-N[H 2 ]\\^ 
The estimates ([T3)l and (fT4|) imply Af : B ^ B and is contractive if 
(15) 



^ + <l±|», e ,„ d! L + ^ 1+s)l|Jfoll<1 



Recalling that p ^ 1, (fT2j) implies that the conditions in (fT5|) are satisfied if e = Thus, (fT0|) has 



a unique solution iJ in a ball of size 2§||i?o||- an d so does the equivalent equation [6]) . Reverting 
the changes of variables, the corresponding y is yt since H £ Si implies the decay of y at ioo that 
uniquely determines yt- 



Lemma 3. With H as in Lemma\^ we have for x 6 Clj, 

1 ,,„„ 1 



H'(x) 



H 



U\x\V 2 9|x| 9 / 2 
Proof. Differentiating ([TO]) and using (fTTj) we obtain 

H(t) , ff 2 (i) 



-HI 



392 



625|a;| 7 /2 



(16) ff'(x) 



cosh(x — i) 



At 2 



2y/t 



dt 



r 

I sinh(a; — t) 

J zoo 



7-392 
2 • 625t 9 / 2 



f/f. 



where the last term is the result of integration by parts. Since ^ t 5 / 2 |jiJ||, | sinh(x — 

t)\, I cosh(a; — t)\ ^ 1, for x,t € Cli in (|16[) . the result follows. 

Remark 2. In order to obtain small error bounds for H and H ' at x = xq := i (24 • 1.7) 5/4 /30 = 
i3.437- • • corresponding to z = z = e t7r / 5 1.7, a good choice is p = \x \ = 3.437- • • and e = ^ 
(for which the conditions in (|15l) in Lemmadhold). With this choice, \H(x )\ ^ |x |" 5/2 ||H|| ^ 



(4 

(40) (3125) 



5/2 ||#o| 



41 784 
40 3125 



~ 5 / 2 . Also, Lemma |31 is applied to bound H'(xq) by using \\H\\ < 



(41)(784) By ©, © and © we have 



2/(2) 



(17) 
Defining 

2/0(2) 

straightforward calculations show that 



25x 2 



; 1- 



T/4 (242) 5 / 4 



30 



25x< 



y(zo) - 2/0(20) = i x j —H(x ), 



■-0 



(18) 
Note also that 

(19) 2/0(20) 



2/(20) - 2/0(20) = i 
ncd bounds 

2/(20) - 2/0(20) 



1 



H(x ) 5 , 

— — + -X H (Xq) 



6 V2oa/^o 

and thus, using the obtained bounds on H(xq) and H'(xq), we get 



890' 



2/' (20) - 2/0(20) 

Cie" 2 "^ 5 , d = 



29 
4468 



-0.5394994- ■• 



(20) 



2/0(20) 



z Q V 12 75a; 2 , 



C 2 e i2V5 , C7 2 = 0.148075 •• 
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Remark 3. The function h is real-valued in Indeed, with x — ir, (O becomes 

392 



(21) 



= o 



By complex conjugation symmetry, if h(r) is a solution of (|21[) then so is h(r); by uniqueness, 
h(r) is real-valued. 

3. The tritronquee in the region x g Ox U fl 2 
Consider the domains 



{x: |x| 



> Po , ^ «S argx 



7T 1 



7T 7T Pq 

— area; ^ — , Re x ^ — = 
4 6 4 V2 



and define the Banach space S2 of analytic functions in the interior of f2i U O2, continuous up to 
the boundary, equipped with the weighted norm 



\H\\ 



sup 

a;eOiUn 2 



,5/2^ 



We consider the operator TV", defined in (|10p. now acting on 6>2- 

Lemma 4. For Po ^ (24 ■ 1.7) 5 ^ 4 = |xo|, there exists a unique solution to the integral equation 
MO]) . H = M [H], in S2, corresponding to the tritronquee through the transformation |17[ ). 

Proof. We first obtain bounds on 
(22) H (x) = ( sinh(x-t)- 

J loo 



392 



625t 7 / 2 



196 

625t 7 / 2 



dt 



196 

625t 7 / 2 



dt 



Consider first x G O2 , The contour in the middle integral in ([22]) can be deformed to a radial one 
joining 00 to x G Oi U O2 for which | e 37 * | $J 1 implying 



(23) 



n 5/2 



» 1 

t 7/a 



2 

< - 

5 



For the last integral in (j2"2")l . no such radial path deformation is possible because of growth of e*. 
On the vertical integration path, we paramerize t = x + ir, t G K + . If we separate out real and 
imaginary parts of a;: x = a + ib, then x G O2 implies -j^= ^ Re x = a and |6| ^ a. Then, 



(24) 



tV 



:dt 



(a 2 + {b + rf) 



2N-7/4 



drr 



1 



,5/2 



{l+p 2 )- 7/4 d P ^ 



5/2 



'-6/0 a 
Combining with (|23|) it follows that for x G O2 we have 

196 



{i+p 2 y 7/i d P 



625 



2 5/4 



{l+p 2 )- 7 ' A dp + 



32 
25 



For x G fii, we note that \t\ ^ |x|, and on a vertical contour of integration | cosh(x — t)\ ^ 1 and 
\dt\ \/2^|i| and therefore 

784^ 



H X) 



J\x\ 



392 
625|t| 7 / 2 



d\t\ sC 



3125|x| 5 / 2 



This is clearly a smaller bound than the one for x G ^2- Therefore, for any x G fii U O2, we have 

32 

\\H \\ sC — =:M 
II oil 25 

For the nonlinear term, the calculations are similar. For x G f2i, on the vertical path, since 
I sinhfx - t)\ < 1, \H(t)\ < \t\-^ 2 \\H\\ and < V2d\t\, 

sinh (. x - t)^§-dt < ^XLlli/ll 2 



2V* 



9|x| 9 / 2 ' 
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For x € 2 , wc split sinh into two exponentials and break the integral accordingly; in one of the 
integrals the contour can be deformed into a radial path. In the other, we parametrize the vertical 
path as in the estimates of H$. Since the bound in f2i is clearly smaller, this results in 

H 2 (t) 



(25) 
where 

(26) N 
Now consider the linear term 



/ sinh( 

J ioo 



1 

18 



x-t)- 

2Vt 

/oo 
[1- 



■dt 



^Np^ 2 \\H\\ 2 , 



-n/4 j ^ 203 
^138 



H(t) 

sinh(a; — t) — ^-dt 



For x 6 fix, as before, we obtain using | sin(a; — t)\ ^ 1 on a vertical path, \dt\ ^ >/2d|t| and 
\H{t)\<\\H\\\t\s/\ 



r 

/ sinh(x — t) 

J ioo 



H(t) 
At 2 



dt 



V2 



Upo\x\ 5 / 2 



\H\\ 



For x £ similarly writing sinh(a; — t) in exponential form, breaking the integral accordingly 
and separately estimating each term we get 



x Hit) , 
sinh(x-t) — -r-dt 



At 2 



L 



Po\x 



5/2 



\H\\ 



(i + p-r'^dp < i 

o 



Clearly, the bound for x G is larger. We conclude that 



r 

I sinh(a; — t] 

J ioo 



m 

At 2 



dt 



^ L \\H\ 
Pa 



Af : B —} B is contractive if for some e > we have 

(27) Lpo l (l + e) + N P q 2 M(1 + e) 2 ^ e , Lp^ 1 + 2N P „ 2 M{1 + e) < 1 

Both conditions are satisfied for p = |^o|i when e = |, i.e. ball size is §M, and the lemma follows 
from contraction mapping theorem, and the fact that fij is the boundary of fli implies that this 
solution is the same in Lemma [3J -corresponding to the tritronquee. 



4. Analysis of y t for x e fLi 



Definition 5. Let 



4 := 



I a; : \x\ 



x\ p ^ 3 , argx £ 



The results in [T] imply that for large x in s.t. 
asymptotic expansion 



Se~ 



•Jx 



7T 7T 

2'~4 
12 



> m > has the uniform 



where £ = , and each Fj(£) is a rational function with poles at £ = 12. Further, since the 

Stokes constant 5 = = i0. 61804.. (cf. [TU], [TT]) is moderately small, £ is also moderately 

small for 1 G U4 even for p = 3. Therefore, a few terms in the Taylor series of each Fj at £ = 
and a fairly small number of Fj's are expected to yield a good approximation of h(x). 
With this expectation, wc choose an approximate expression ho for h t in SI4 in the form 



(28) hoi x)=\t + ^ + ±- + i- + 



e 

6 



i 3 
48 



432 20736 / 



11 
72 



43 

1152 



9i 



128a; 2 
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where 

< 29 » «=^ s ='^ 

Let 

(30) h(x) = ho(x) +x- 1 ' 2 G(x) 
By (gl), G{x) satisfies 

(31) G" (1 + ho(x)) G(x) = _ R{x) ^ 
where R(x) is given by 



R(x) = ^ U» + h' - ho(x) - \hl{x) - Jgj) = £ x-V* rj (e x ) 

V ' 3=5 

and the rj(Q are polynomials in C _1 i where only r-j has a nonzero constant term — 1|| (see 
(|143p ~ (|146p in the Appendix for the precise expressions of tj). Define 



9 

(33) Rq(x) = x- 5 / 2 r 5 (e x ) , R^x) = x" 3 r 6 (e x ) , R(x) = Y^x~ j/ ' 1 r j {e x ) =R-R Q -R 1 

3=7 

and 

392 

(34) f 7 (C) = r 7 (C) + ' ^ (0 = ^ (0 for i + 7 



To write (|31() in integral form, we need the properties of the Green's function of the operator on 
the left side of the equation. It is more convenient to write a nearby equation with an explicit 
Green's function, and for this end we find quasi-solutions of the homogeneous equation 

(35) u" - (1 + h (x)) u = 

Formal asymptotic arguments for large x suggest that one solution of (|35l) has the asymptotic 
behavior 

J(x) 



(36) u ~ yi (x) = e~ x I 1 + 

where 

_ Ser x S 2 e~ 2x \§Se- x SVj^ 5S 2 e~ 2x 25S i er ix 

( ' ^ X >~ 3 + 16y/x 72a; + 108x 48a; 3 / 2 + 20736a; 3 / 2 

We can readily check that y\ solves (|36p up to 0(x -5 / 2 ) errors: 

(38) 2/1 -(l + fco(i))m =q(x)yi, 
where 

9 

(39) <z(*) yi (x) = 5>-^(^), 

3=5 

Here, all qj(Q are polynomials in l/£ of degree at least 2 (see equations (|147|) - (|150[) in the Ap- 
pendix). We chose j/i to ensure that the error term q\y\ is 0(x~ 5 ^ 2 ) for large x G O4. 

A second independent solution to the homogeneous equation (|3"5|) 2/2 is given by 2/1 (x) j J a ^^tj | • 
With a suitable choice of integration constant, 2/2 becomes 



2 



^\og(ix) + z 2 (x) 



(40) 2/2 Or) = ift(a:) 
Here 

e 2x 2Se x 

(41) 22,0(2;) = — . z 2 ,i{x) = ~t^J=- 



, where z 2 (x) = z 2 , {x) + z 2 ,i(x) + z 2 , R (x), 
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(42) 



dx 1 



2St 



Se x 



5S 2 



3VX 7 3a;' 3 / 2 24x' 



Using the fact that y\ and y 2 defined above solve (|38[) and their Wronskian is y\y' 2 — J/2 y'\ = U 
inversion of pip results in the integral equation 



(43) G(x) = G (x)- 

where 
(44) 
where 
(45) 

(46) 



{V2{x)yi{x') - y x {x)y 2 {x')} -V(x')G(x') 

V(x) = V (x)+q(x), 
V (x) 



G 2 (x'p 



dx' =: Af[G] (x) 



1 

Ax 2 



Go(x) 



{yi(x)y 2 {x) - y 2 (x)y 1 (x')} R(x')dx' 



Definition 6. Let S4 be the Banach space of analytic functions in the interior of ^4, continuous 
on ^4, equipped with the norm 



(47) 



||G|| = sup x 5/2 G(: 



The usual sup norm will be denoted by 



We seek a solution to (|4"3"j) . i.e. G = N [G] in S 4 with p > 3. It will be proved that N : <S 4 — s- 5 4 , 
see ((33]) . The general integral reformulation of (j3"Tj) 



(48) 



G(x) = C m (x) + C 2 y 2 (x) +N[G] (x) 



implies (Ci,C 2 ) = (0,0) since neither yi nor y 2 are in £4. Therefore, any solution to (f3"Tj) in 1S4 
must necessarily satisfy P5|) . 

We now prove the following result: 

Theorem 2. For p ^ 3, £/iere exists unique solution G to the integral equation J^3[ ) in a ball 
-B4 C ^4 of radius 4- Through the change of variables in ([3]) and (|30[) , </i?s corresponds to the 
tritronquee solution y = yt- 

The proof of Theorem [2] follows from the following Lemmas that are proved in subsections 

USES 



Lemma 7. for p ^ 3 we have 
(49) 

Lemma 8. For p ^ 3 ; and G £ £4, 

(50) 



IIGoll <2 

(y2(x)y 1 (x / ) - yi(x)y 2 (x')) V(x')G(x')dx 
Lemma 9. For p ^ 3, and G, G\,G 2 E S4, 
(51) 



1 



{V2{x)yi{x) - yi{x)y 2 {x')) —^=G 2 {x')dx' 

2vx' 



^IIGII 2 , 



(52) 



(y 2 (x) yi (x') - yi (x)y 2 (x')) —= [G\{x') - G 2 (x')] dx' 

2\/x' 



^(IIGrl 



|G 2 ||)||Gi-G 2 | 
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Proof of Theorem [2] It is clear from Lemmas M and M that M : B 4 ->• B A (cf. ([35])): for G, 

Gi and G2 in S4 we have 



and 



\W[G] j| ^ ||Go|| + \\\G\\ + ±\\Gf < 2 + I + g < 4. 
HA^Gx] - TV [G 2 ] || < |||Ci - G a || + A|| Gl _ G 2 || < ^||G a - G 2 | 



By the contraction mapping theorem, there is a unique solution to (|43p in B4 if p ^ 3. From 
([30)) . it is clear that G corresponds to a solution /i of (|4]) that is singularity- free in the closed 
domain fi 4 and has the leading order asymptotic behavior h ~ s ^ as x — > —ioo on the negative 
imaginary axis. By [T] §5.2 (see also [3]), for any G there is a unique solution h oi ^ with the 
behavior Ce~ x x~ 1 / 2 as x — > —ioo analytic for large i in a sector in the fourth quadrant |W| The 
value G = £ identifies this solution with the tritronquee. (This also follows from classical results 

cf. mo 

4.1. Preliminary Lemmas. In this subsection we obtain various integral estimates needed in 
the sequel. 

Definition 10. Let P be a polynomial, P(rf) = J2m=m PmV~ m - We define the following weighted 



norms: 



2 2 

\Pm\ for j > 2; J" 2j [P] = ^ —\Pm\ for m > 



.7-2 



9 mi mi i 2 + 2-j - 2 

^•[P] = ^ E Wfor j>3; J- 4J [P]= •( U 

7 — 3 z — ' 7(7 — l)m 

m—rno m— mo 7 

Lemma 11. //Zo > 2, g is analytic in ^4 ™£/i ||<?||oo < OO; and 

L 

(53) u,(aO = J2 x ~ W2p i ( eX ) - 



bm| for j > 1 , m > 



Z = Zn 



where m) = E™L *WC~ m the 



(54) 



where 



[ X g(x')w(x')dx' <|MUX>l"' /2+1 Q'> 

./ —zoo , , 

i— in 



-//2 



and thus the 



Proof. The various terms in the integrand in (|53|) are of the form g(x')p mi ie~ m: 
contour of integration can be deformed to a radial path from 00 to ir € Sl 4 and that |e _ma: <7(;r)| ^ 
1 1 .g 1 1 oo for a; e il. Then, clearly, 



Pm,l / 9(x')e 



- mx 'x'~ l/2 dx 



< bm.illlsl 



.1" 



n-i/2 



d|x'| < 



|p m ,z||x|-'/ 2 + 1 
Z/2 - 1 



I 



Lemma 12. If Iq > and 

(55) 



where Pi(() = ^2m'=iPm.iC m * s a polynomial in £ , then 



(56) 



/ u>(:r/)dr/| ^ 
"'- 100 z=z 



-1/2 



analysis in [l] is done in the first quadrant, but by symmetry w.r.t. 1 £ 1 it applies with straightforward 
modifications to the fourth quadrant. 
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where Qi = F 2 ,i [Pi] 

Proof. We note that the terms in the integral in (|56|) are of the form 

,-l/2\ lp m ,l _ 



p m ,ie~ mx V 1/2 



dx' 



Pm,l c -mx' x 

m 



2m 



3 rax I l /' 2 1 



Therefore, integrating out the first term explicitly and deforming the path for the second term to 
a radial contour, it follows that 



Prn,l 



e- mx 'x'~ l/2 dx' 



2\prnJ 1/2 



m 



I 



Lemma 13. If for l > 3 w(x) = Y,Li x l/2p l( eX )> where Pi(() = EmLoZWC m > and */ 9 



analytic in with \\g\\oo < 00, then 

f log - g(x')w(x')dx' < \\9\UJ2Q 



-1/2+1 



l=l 



where Qi = J" 3 ,/ [Pi] 

Proof. Once again because of the analyticity and decay of integrand, we may deform the inte- 
gration path to a radial one joining 00 to a; £ O4. The general term in the integrand is of the 
form 



Pm ,ix'- b/ e- mx g{x')\og-dx' 

x 



f 1/2 ^^ mX ' rJ 

Since it is readily checked that for \x'\ 5^ [x 

X 



log — 

x 



1/2 



it follows that 
(57) 



p m .ix'- l/2 e- mx 'g(x')\og-dx' 

x 



\oo\Pm,l\\ x \ 



-1/2 



■r 



n-Z/2+1/2 j| _/ 



d\x'\ = 



?Wi 1-1/2+1 



1-3 



I 



Lemma 14. ///or Z > 1 ; = £f =io x~ l ^ 2 Pi(e x ), where P^Q = Y2=iPm,lC m , then 



log — w(x')dx' 

x 



Z = Z 



-J/2 



w/iere = J4,z 

Proof. Once again because of the analyticity and decay of the integrand, we may deform the 
integration path to a radial one joining 00 to x <= Q4. A general term in the integrand is of the 
form 



Pm,ix'- [/2 e- mx ' 



log- 



dx' 
I 



Prn.l 



in 



J/2, X \ , Pm,l 



x'-'^lo 



X 



-1/2-1 



■ p m!ie - m -'x'- l/2 - 1 \og- 
2m x 



Noting that the complete derivative is zero at the end point x' = x and applying Lemma [13] to 
bound the integral of the third term (on a radial path) we immediately obtain 

'2 I 



> x' 1 

Pm,ie mX log— Cfo' < —\p m ,l\ 

x m 



I l-l 



-1/2 



from which the Lemma follows. 
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4.2. Bounds on J, yi, Zi r and z 2 for x € CI4. Let 

S 



(58) 
and define 



yi,o(x) = e x , 2/1,1 (x) = — -=e x , y 1M = 2/1 - 2/1,0 - 2/1,1, 



35 



(59) = 
comparing with (|37|) we see that 
(60) 

Note that for x G O4, 
(61) 



19 S 2 e- 2x 



30 



..-1/2 



_5_ „ - x J5_ 
"16 + 6912 



5 3 e 



J(x) = ^ (l ■ 



, / \ 1 3151 19 1 |5| 2 5 151 25I5I 3 



10 24 Jp 36^p 16 p 6912p 
Using J (a;) in (1571) . it follows that for x € SI4 we have 



(62) 

From (|36p. it follows that 
(63) 

Now, {37]) and (J5HH imply 
(64) 



e^x) 



=: Jm 



< 1 



M 



xe 2x yi^ R (x) 



From (|58j) . it also follows that for 3; 6 ^4, 



(65) 



e x (2/1,0 (x) +yi,i(x)) 



^ 1 



Expressing j/i in terms of J, as in (ptBf and (J37J, in (|4"2"j). it follows that 



(66) Z2,r(x) = Z2,r,o(x) + / E(x')dx' + — 



36 / r' 3 / 2 

OVJ 1/ — zoo X 



J(x') 



where 
(67) 



^2,_R,o(a;) 



235 2 3615 2 235 3 e 



3„ — x 



-->77S 4 e~ 2x 



72x 3456x 2 216x 3 / 2 41472x 2 ' 



3„-3x \ x -x 



2J(x') 3J 2 (x') 



dx', 



(68) 



E(x) = ^x-^E j (e x ), 
3=5 



where each £j (£) is a polynomial in l/£ with no constant term (the precise expressions are given 
in (|151[) - (|152l) in the appendix). We note that 



(69) 



XZ 2 ,R,0 



^ 23|5| 2 + 23|5| 3 + p61|5| 2 + 577|5| 4 j 



72 216pV2 ^ 3450 41472 



Using Lemma [T^J it follows that 



(70) 



/x 8 
E{x')dx'\ < •' : /,„, =: E M 

' ica j=5 
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(see Definition ITUl and the expression of Em in the Appendix, (|153p ). On integration by parts, we 
get 



(71) 



7Se~ 



3(> 



,3/2 



dx' 



7S 7S 



36a; 3 / 2 24 



P /5/2 



Therefore for x £ we get 



(72) 



7Sxe~ 



3(i 



r /3/2 

ZOO ^ 



da;' 



7|5|(>/2 + l) 



where we used the fact that on a vertical contour connecting — ioo to a; £ Q4,, \dx'\ ^ \/2d|a;'|. We 
note that 
(73) 



x \ l + x-^j) -l + 2x- 1 / 2 J-3x- 1 J 2 [ =-4e 2x x- 3 / 2 J 3 + 



hx~'e ix r 



-5/2 e 2xj5 



4S 3 e 



4S 3 



and 

(74) -4e 2 "a;- 3 / 2 J 3 (a;) - 

y ' y 1 27a; 3 / 2 27a; 3 / 2 

Now, using Lemma [12l it follows that 



l + x-V 2 J {l+x- l / 2 Jf 



(75) 



x AS 3 e 



-dx' 



o 27a;' 3/2 

Deforming the contour to a radial one, it is clear that 



8|S| 3 
27pV2^ 



(76) 



45 3 x 
27 



?'oo X 



,3/2 



1 



i(z') 



1 



da;' 



4|5| a 
27 



3j, 



? 2 

VP p 



Therefore, using ([75]) and ([76]) in ([74]) we get 



(77) 



4e 2 *'j 3 (a;') , , 8|S| 3 4|S*| 3 J 



,3/2 



-dx' 



-IQG X 

From ([73]) and ([77]), it follows that 



27 VP 



Jm Jm_ 

Vp 3 p 



(78) 



2x ' (l + x'- 1/2 J(x')^ 2 - 1 + 2x'- 1/2 J(x') - 3a;'" 1 J 2 (a;')dx' 



8|^| 3 4|S| a j, 
27VP 



9 



1 



Jm J r 



0J M 



2p- 1/2 Jl 



M 



\/~P ZpJ 1-P~ 1/2 Jm 3(l-p- 1 / 2 J M ) 



Therefore, combining all the estimates ([59"]) . ([70]) , ([72"]) and ([75]) in the expression of z 2 ,fl in ([BTi]) 
we get 



(79) 



23|Sf 23|Sf /361|Sf 577151' 



72 



216/9 1 / 2 V 3456 41472 



p- 1 +E M 



7|5|(V2 + 1) | 8|S| 3 ^\S\ 3 j m ( 1+ Jm + Jl\ 5Jj, | 2p-V*J* M 



36VP 27VP 9 
Therefore, using ([4*0]) and (f4"T]) we get 

(80) 



3^; i-p-Wm 3 (1 _ p -i/2 Jm ) 



2 — : Z 2,R,M 



-2x I ^ 1 , I -x 1 , Z 2,R,M 1 . 2|5| Z 2,R,M 



P 



2 3VP 



/0 



To help the reader who would like to check the intermediate steps in the calculations, we provide 
in the Appendix the numerical values for p = 3 of the various constants appearing in the estimates. 
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4.3. Bounds on V(x) and proof of Lemma [8J We first seek bounds on q yi . It is clear from 
the expression of qy\ in (|39f that Lemma ITTI applies to w(x) — x~ 5 ^ 2 e 2x qyi, g(x) — x 5 / 2 G(x); 
noting that ||G| = j|x 5 / 2 G|| 00 , we obtain 



(81) 



14 



x 5 / 2 / e 2 *'q(x>) yi (x')G(x')dx> < \\G\\ £ 7^ fe_ 5 ] p^/ 2 + 7 / 2 = M JG|| 

J-ioo j = 1Q 



(see Def. [T0|) . The explicit formula of M q is given in (|1 54|) in the Appendix. Further, 
(82) x 5 / 2 yi (x) [ (z 2 (x')- z 2 (x))q(x') yi (x')G(x')dx' ^2Y ltM Z2,MM q \\G\\ 



In ([82]). recalling that q(x)yi(x) is a polynomial in l/e x of degree at least 2, we applied Lemma 
[IT1 with g(x') = e~ 2x ' z 2 (x')x /5 ^ 2 G(x / ) and w(x') = x'~ 5 ^ 2 e 2x> ' q(x')yi(x') in the part of the in- 
tegral involving z%{x'), while in the second one, we took g(x') = e x ~ x x' 5 ^ 2 G(x'), w(x') = 
x' 5 / 2 e x qi(x')yi(x') and used e~ x z 2 (x)y\(x) € z-> atYi m. Lcmma[T3l for w(x) = x~ 5 / 2 a(x)vi(x), 
g(x) = x 5 / 2 G, implies 



(83) 



5^ 
24 



r 5/2 



yi (x) f \og X - q(x') yi (x')G(x')dx' < ^f\\G\\ £ F 3tj [q^p-i' 2 ^' 2 

J— ioo X ^ :_in 

. 5|5| 2 



j=io 



24 



M L JG\\ 



where the detailed expression of Ml, q is given in the Appendix, (|157|l ). Therefore, using (|4"0|) , 
(I82l) and (El), it follows that 



(84) 



(2/2(x)yi(x') - yi(x)y 2 (x')) q(x')G 



2z 2M M q + ^~M L>q 



We now bound the terms involving Vq — jp-. Noting Lemma [Til applies to w(x) = e x Vqx b l 2 
and g(x) = x b ^ 2 e x Gy\\ since both ||G|| = II^^GHoo and \e x yi \ < Y^ M , it follows that 



(85) 



a- 5 / 2 / V r (x / )»i(a/)G(a/)d!B' 



Since e^j/i, e 2a: 2!2 and x 5 / 2 G are bounded by Yi,m, ^2,m and ||G|| respectively, we have 



(86) 



x b ' 2 e- x 



V {x') yi (x')z 2 {x')G{x') 
^Yx,mZ2,m\\G\\ (m 5/2 



-x-\-x I 



U\x 



where we used the fact that on a vertical contour joining —ioo to x € f^4 we have \e x x \ = 1 and 
| da/ 1 < V2d|x'|. Lemma li"3l applied to iu(a;) = e~ x a;~ 5 / 2 Vb, g(x) = e x yiX 5 ^ 2 G gives 



(87) 



J — ioo x 



Therefore, combining (|85]) -(|87 p and using (j40|) . we obtain 



x 3/2 / (l&(a0tfi(aO - l/ifr)^)) V (x')G(x') 



A/ 



, 5|5| 2 
^4T Z2 - M + _ 



l|G|| 
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Collecting the contributions of the terms involving q and Vo in (|84[) and (|55|) respectively, it follows 
that 



(89) 



(y 2 (x) yi (x') - yi(x)y 2 (x')) V{x')G{x')dx' 

O 

Y X>M I [2z 2 , M M q + ^-M L ,^j + Y hM 



y/2 + 1 5|S| 2 

-Z2,M + 



14/5 



288/? 



l|G|| 



Vm||G|| 



Since all the quantities involved in Vm are decreasing in p, it is clear that for p ^ 3, Vm is bounded 
by its value at p = 3, which in turn is less than 9/40 and Lemma [8] follows. 

4.4. Nonlinear terms and proof of Lemma [9J Applying Lemma [TT1 to w(x') = e~ x x' -11 ^ 2 , 
g{x') = ^e x ~ x x' 5 G 2 (x')e x yi(x'), noting that je 31-1 ^ 1 on a radial contour in O4, and finally 



that |x' 5 G 2 (x')| ^ ||G|| 2 , \e x ' yi (x')\ Fi M , we obtain 



2 l„ai',,. 



(90) 



e x j/i(x)e- 2x z 2 (x)x 5/2 



_ x ,e x ' Vl (x')G 2 {x') ^ , ^ F 2 M z 2>Af 



2^ 



-dx 



l|G|| 



Furthermore, we note that 
(91) e x yi (x)x 5 ^ 



. xW e x ' yi {x')e- 2x ' z 2 {x')G\x') ^ 
2VX 1 



9p 2 

V2Y* M z 2 ,M 
9p 2 



\G\\ 



where in (|9~Tj) we have |dx'| $S v2d|x'| and \e x x \ = 1 on the vertical contour joining —ioo to 
x £ Q4. Applying Lemma [TBI to w(x) = e~ x x~ 11 / 2 , g(x) = x 5 G 2 e x yi(x), we get 



(92) 



5tP 
24 



yi(x)x 



5/2 



log f= ax 



x 2VX 7 



5\S\ 2 Yl 



M M/OM 2 



Combining (|90|t. (|9T|) and ([92]) we obtain 
(93) 



(y 2 (x)yi(z') - 2/i(x)y 2 (x')) ^©dx' 



2Vr 



1 l.M 



\m P 2 



1 + V2 
9 



l|G|| 



^2.Af 



192 



J l|G|| 2 
7m||G|| 2 



A very similar calculation shows that 

(y 2 (x) yi (x') y,(x)y 2 (x')) 9l^l^Ml dx ' 



<T M (||Gi|| + ||G 2 ||)||Gi-G 2 | 



Proof of Lemma [9] follows since a calculation of Tm, which is a decreasing function of p, 
18 ^ j_ 

467 v 25 



shows T A/ < ^= < ± for p ^ 3 



4.5. Bounds involving Go and proof of Lemma Using (f3"2" j) - (f3"3")l and the form of y 2 in 
(l4"0l) . it is convenient to decompose Go defined in (|4"6"|) as follows 

(94) Go = Gq,i + Go, 2 + Go. 3 + Go,4 + Go, 5 + Go, 6 + Go, 7 
where 

(95) G ,i(x) =yi(x) ( y x {x') [z 2 {x') - z 2 {x)] R(x')dx' 

J — ioo 

(96) G 0>2 (x) = yi (x) [ y x , R {x') [z 2 (x') - z 2 {x)] (R Q (x') + R^x')) dx' 

J — ioo 

(97) G , 3 (x) = yi (x) [ [yi,o(x') + y M (x')] [z 2 , R (x') - z 2 , R (x)} (Ro(x') + R^x')) dx' 
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(98) G 0>4 (x) =yi(x) / [yi,o(x') + yi,i(x')] {z 2>0 (x') + z 2 ,i(x') 

o 

-z 2fi (x) - z 2 ,i{x)} (R a (x r ) + Ri{x')) dx' 



(99) 
(100) 



5S 2 

Gq j5 {x) = -t— yi(x) / log — yi(x')R(x')dx 



24 



592 

Go,a(a;) - -24-Wi(!c) / log - Vi,r&) + ^lOOl 



KC2 /-x / 

(101) G 0t7 (x) = -=- Vl (x) log-(y 1 , (x') + y 1 , 1 (x , ))(Ro(x')+R 1 (x , ))dx' 

^ J —too x 

In £)4.5.1l below we obtain bounds Mj for ||Goj|| for j = 1, 2. .7; using those, we get 



(102) 



\G \\^J2 M i 
j'=i 



The formulas for Mj,j = 1, 7 are given in the following subsections. These expressions will be 
shown to be decreasing in p, and Lemma [7] will follow using the values of Mj at p = 3. 

4.5.1. Bounds on Go,i- Using Lemma [Til with g(x') = e x ~ x ' e x ' yi(x') and w(x') = e~ x ' R{x') we 
get 



(103) 



p x 9 

.T 5 / 2 yi (x)z 2 (x) / yi (x')R(x')dx' < FrV^.M X] ^ fa] P~ j/2+7/2 

J— ioo j—7 



Note in the estimate above, the factor e x outside the integral was placed inside the integral in 
g(x'); this is legitimate since in a radial deformation of the integration path from oo to x £ we 



Noting again that R{x) + 62 ^t.y 2 has only decaying exponentials, Lemma [IT] (this time with 



g(x) — e x y\(x)z 2 (x) and w(x) = e x (r(x) + |||a; 7 / 2 ^ implies 

392 



(104) x 5/2 yi (x) / y x {x')z 2 {x') 



R(x') 



625x' 7/2 



dx 



^Y 2 M z 2 ,u^iAfj]p- j/2+7 

3=7 



/2 



Since on a vertical contour in joining —ioo to x G SI4, |e x+x \ = 1 and \dx'\ ^ \/2d|£'|, we get 



(105) 

It follows that 



x 5 / 2 yi (x) / yi(x')z 2 (x') 



392 



625x' 7/2 



784 xr2 r- 
^ 3125 Y ^ mZ2 - mV2 



(106) ||G 0l i|| < Y? M Z2, M { ^^2 + ^-7 [rr + h] + 2^ Tij [r,] }> =: ^V^.m^i =: M x 



where the explicit expression of Mq,i is given in (|164[) in the Appendix. 

4.5.2. Bounds on Go. 2- From (|64j) . f|80[) we note that e 2x xyi,R and e _2:c 2:2(2:) are bounded by 
^i,_r,m and 22, M respectively. Lemma [TT1 applied to w(x) = A (i?o + i?i) implies 

8 

(107) ||G , 2 || < 2Y 1 , Af z 2 , M Y 1 ,^ Af ^ p-J/2+V2 3 [ rj ._ 2 ] =: 2Fi, Af 2 ; 2,Mli 1 R,MM Gj 2 := M 2 
where the formula for M(j i2 is given in (|168[) in the Appendix. 
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4.5.3. Bounds on Go, 3' Note that xe~ x Z2,r and e x (j/1,0 + Vi,i) ar e bounded by zi,r,m and l+<j-^= 
rcsp. Lemma [TT] applied to w(x') = [Ro(x') + R\(x')] jx' for the term containing Z2,r{x') and to 
w(x') = [Rq(x') + Ri(x')] for the one containing Z2.r(x) implies 



(108) 

where 

(109) 



l|Go, 3 || <Y ltU z 2iRtM 1 + 



3VP 



G,3 



A/3 



M G ,3 = M G , 2 • E'' ' fo] . 

j=5 



The concrete expression of Mq^ is given in (|169p in the Appendix. 
4.5.4. Bounds on Gq^: Using (|3"3"|) and (|5"5|) . it follows that 



(110) 



T(z) =: (j/ li0 + 2/1,1) (i?o + = E *~ i/a *i(e X ). 

J=5 



where tj{C) are polynomials in 1/C, having no constant and linear terms; the precise expressions 
are in ([r7D|) - p7T|) in the Appendix. 



(Ill) 



U{x) =: (yi,o + yi,i) {z2fi + Z2,i) (R0 + R1) = ^2x~ j ^u J (e x ) 



3=5 



where Uj(C) are polynomials in 1/C without constant terms-see (|172p ~ (|173p in the Appendix. We 
also note that 



(112) 



T{x) 



dx 



E x ~ j/2 r i ^) 



3=5 



+E 



3=7 



where Tj(C), tj(£), are polynomials in 1/C, see (|174p - (|177[) in the Appendix, with no constant or 
linear terms. Again, we note that 



(113) 



U(x) = — 
dx 



3=5 



10 



E* ' L V, ;«•••; 

3=7 



where the polynomials in 1/C, i*j(C)i v j{0 have no constant terms, sec (|178|) - f)181[) . Using (|98l) . 
PTu| - (fTTgi) . it follows that 



(114) G 0A (x) = Vl {x) I J2^ j/2 ^3 (? X ) - Mx) + z 2A (x))J2^ 3/2 r 3 (e*) 



3=5 



3=5 



r x { 10 9 

J-ioo [ J=7 j=7 



Straightforward calculations show that 

(115) E x ~ j/2 »i ~ Mx) + z 2 ,i(x)) E x~ 3/2 r3 (e x ) = £ ^'^3 (?) , 



3=5 



3=5 



3=5 



dx' 



where Pj(C) s are a l so polynomials in 1/C with no constant term; see the Appendix, starting with 

eq. CEH2D. 

Applying Lemma [TT1 to the two terms in second integral on the right of (|114p . using \e x ~ x | ^ 1, 
w(x) = Y^j=j x ~^ 2 tj (e x ) for the integral on the first line and w(x) — ^7=7 x~^ 2 Uj (e x ) in the 
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second line we get 

(116) ||G ,4|K^l,M 



5> 

3=5 



3/2+5/2- 



10 



+ ^Zp~ j/7+7/2 J r i,j [%] }• =: *i,m (A/gao + M G)4 ,i) =: M. 

3=7 



*4, 



see ()184|) . (|185|) . where pj is the polynomial obtained from pj by replacing each coefficient by its 
absolute value. 

4.5.5. Bounds on Go. 5, Go, 6 Go, 7. For Go, 5 we simply use Lemma [T3l with w(x) = e~ x R(x), 
and g(x) = e x y\. We obtain 



(117) 



5/2 



log- yi (x')^(^) 



24 



J =7 



24 



where Mg^ is given in (| 186[) in the Appendix. We again use LemmafTBIwith w(x) — ^— (Rq + Ri) 
and g(x) = e x xyi^u, to obtain 



(118) x^ 2 ^- yi I log-tfi lfl (a/)(flo(a/)+iZi(i / )) 



5S 2 



24 



<5|5| 2 A 

^ -^:7-'l,AfJl,Ji,Af 2 



.-372+7/2 



•^3,j [rj- 2 ] 



3=7 



5J5P 
24 



ii,Myi,fl,MM G ,6 =: M 6 , 



where Mg,6 is given in (|187|l in the Appendix. Now consider Go, 7 (x). Recall that 

7 

(119) T(x) = (y 1>0 + (jRq + fix) = 5^ aT^ (e») 

3=5 

From the expressions of tj in (|170[) - (|17ip . it is clear that (|119[) involves only decaying exponentials. 
Therefore, we may apply Lemma [T4l to w(x) = T(x) giving 

(120) x 5 / 2 ^- yi (x) \og-{y 1 , Q {x')+y 1 , l {x')){R < >{x l )+R 1 {x'))dx l 

^ J —too x 

s\s\ 2 



24 



-y 1m J2p- 3/2+5/2 



3=5 



5IC-I2 



where Mg,7 is given in (|188l) in the Appendix. 

5. Analysis of y t for |z| < 1.7 

For the analysis of the inner region \z\ < 17/10 it is convenient to use the symmetry of the 
equation and write it as 

(121) g" = 6g 2 + t where g(t) = e 27Tl/5 y{-te™ /5 ) 

We take initial conditions close to (fT9"|) . (|2"0"]) , converted into conditions for g(t) 

(122) g(t ) 



280 150 

'519 ; 3(<o)= 1013 ; wherei o = -L7 



We first construct a polynomial which is sufficiently close to g in L°°([£o,0]) so as to be able 
to rigorously preserve error bounds of the order of those in (|18l) . The way the polynomial go 
below was obtained is essentially by projecting g, calculated from its truncated Taylor series, on 
Chcbyshev polynomials of order seven, enough for the aforementioned goal; the polynomial is 

, . . . 280 150s 239s 2 110s 3 32s 4 9s 5 16s 6 8s 7 

(123) q (t) = 1 1 1 1 1 , where s = t- t . 

\ ) yuw 5ig 1013 10331 1477g Q853 43g7 3Q505 49105 ' u 
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Definition 15. For a polynomial P{x) = X)l-=o c k xk on an interval I we define an I 1 norm by 
ll-Plli = Efe=o \ c -k\m k where m = supj \x\. 

Taking g = go + S we get 
(124) 5" - I2g 5 = 6S 2 + R; 5(0) = o x = .g(0) - 5o (0), 5'(0) = a 2 = g'(0) - g' Q (0) 
where —R = g' Q ' — 6g^ — t is a polynomial of degree 14. 
Proposition 16. For t G [i ,0] we ftcroe |i?| < 1/8619. 

Note 4. fa,) Estimating rigorously and with good accuracy, a polynomial P(x) on an interval 
I = [a, 6] is elementary, and it can be done rigorously and efficiently in a number of ways. The one 
used here is the same as in [4|. We choose a suitable partition of [a, b], II = (xq,xi, ...,x n —\,x n ), 
where xq = a, x n = b; then write t = h(xt + + u on each subinterval [xj_i,Xj] for i = 1, ..n 

and re-expand P to obtain a polynomial in u. The polynomial in u is estimated by taking the 
cxtremum of the modulus of the cubic subpolynomial and placing || ■ ||i on the rest. In practice a 
small number of partition points suffices to get a good bound. 

Since all partition points are nonpositive, to simplify the writing we denote by —II the partition 

(-Xo, -Xl, ~X„). 

(b) Bounding rational functions with real coefficients reduces to estimating polynomials since 
the inequality \P/Q\ < e with Q > is equivalent to the pair of inequalities P — eQ < and 
P + eQ > 0. However, in our case the denominator is W which is very close to one, so an upper 
bound of the modulus of the numerator and a lower bound of W give an equivalent accuracy. 

Proof. Use NoteH (a) and the partition n = fj, §, ^, ^,0). | 

To write the equation for S = g — go in an integral form suitable for a contraction argument, 
we would need the fundamental solution of the linear operator on the left side of (|124[) . Of course 
this cannot be done in closed form; once more, we use a pair of polynomials close enough in L°° 
to a fundamental system and estimate the errors introduced in this way. The pair of polynomials 
is obtained in the same way as go was found, and is given by 

,,„„, T/ , 9489s 2 1350s 3 359s 4 1526s 5 708s 6 503s 7 211s 8 , 

125) Jj (t) = 1 1 1 1 (s = t-t ) 

v ; w 2932 4721 199 3719 1633 2201 6486 v ' 

48s 2 2941s 3 675s 4 1832s 5 2305s 6 677s 7 1573s 8 531s 9 . 

J 2 (t) = s 1 1 1 (s = t-t ) 

w 659797 2730 4873 4745 19401 14054 53783 128216 v ; 

Let W — J\ J 2 — J%J'x be their Wronskian. The equation that the pair (Ji, J2) satisfies is 

(126) /" + Af + Bf = 
where 

(127) A -= w ' B - = w 

The Green's function associated to (|126[) is 

(128) Q(s,t) - M^(s)- 1 [J 2 (t)J 1 (s) - J 1 (t)J 2 (s)\ 
We define the linear operators K\, K 2 such that 

(129) K\ [/] (t) = j* g(s,t)f(s)ds, K 2 [f](t) = jT ^Hf( s )ds, 

We can now rewrite now ()124[) as an equivalent integral system. 

Let B 1 = 12. 9 o + B, S = (S, 6') and r(S(s), s) = AS' + Brf + 6S 2 , J 1 = (J U J{), J 2 = ( J 2 , J' 2 ) 
and vector operator K defined by K[/] = (Ki[f],K2[f]). We have 

(130) S = aiJi + a 2 J 2 - Ki? + Kr = NS 
and (|18| implies 

(131) [ai| < ai := 1/290 and |a 2 | < a 2 := 1/152 

To analyze the integral system we first estimate the various quantities in (|130|) . 
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Proposition 17. The following estimates hold in the sup norm on [to, 0] 

(132) m^iW.hlW.h/Wltt^p max{||J 2 ||,||J 2 /^||K^ ||J(||< §, II4II < § 

(133) |W - 1| < 1/500, ||A|| < 1/1216, ||Bi|| < 1/492; 

(134) sup H01J1 + 02J2II < 1/180; sup ||ai J[ + a 2 J' 2 \\ < 1/90 

|ai|<ai,|a2|<a2 |ai|<ai,|a 2 |<a2 

Proof. The proofs are based on Note [4] (a) to estimate the polynomials and rational functions. 

We illustrate the calculation on a rational function, B\ on a sample interval, say (— |,— yg), 
see below. We thus take s = — 13/10 + u/5 and simplify the resulting expression. Denoting by Ej 
polynomials with i 1 norm less than 1/1000 we simply get on this interval 

/ 3 2u 22u 2 6u 3 \ 

B x = + + E x / 1 + E 2 ) 

V 2 332 60137 3361 7241 J /v ; 

The derivative of the cubic has one root for u = —[1,1]; the value of the absolute value cubic 
there is < 1.3 ■ 10~ 2 . Checking it at the endpoints of the interval as well, we see that this is its 
maximum absolute value. The other calculations are as straightforward as this, so we naturally 
omit the details. 

We found it convenient to use a different partition n for estimating maximal absolute values 
of each function. We chose partitions — |,0) and — (jg, yfj, yg, 0) for J\ and J[ respectively 
(cf. again NotcS (b)); -(g, ±±, i,0), -<g, 1±,0) and -<*g, g,0) for J 2 , J' 2 and W respectively. 
Finally, we use Notc|H(b) to estimate A and B\ using the partitions — (t5' t§' ^' ^R)' ^3' TTJ'^) and 

-(to"'§>TO> TO>°) respectively. 

For (|134p we note that /(s, a) := ai Ji(s)+a 2 J 2 (s) and /'(s; a) are linear in (01, a 2 ) and thus for 
any s, the extrema of \ f{s; a)|, \ f'(s, a)| are reached on the boundary. Thus we only need bounds 
when a\ = ot\, a 2 = ±a 2 . The partition points chosen are = — (yjj, ^,0) for f(s 7 a), — (^X, |,0) 
for /(s; 0:1,-0:2), -(to'TO' ) for /'( s 5 a ) and _ (t5' I> I5>°) for /'( s ' a i' _a 2))- I 

To analyze (|130[) we use the norm |<$||^' = max{||5||, |||<5'||} where ||/|| is the sup norm on 
[*o,0]. 

Proposition 18. (i) The nonlinear operator N in (|130|) is contractive in the disc {5 : ||<5||^ 
1/158}; the contr activity factor is < 1/6. 

(ii) We have e x = \\5 - ai J x - a 2 J 2 \\ < 1/1500. e 2 := - a x J[ - a 2 J' 2 \\ < 1/658. 

(Hi) Let a = 87/469 and b = 41/134. We have 

(135) \g(0) + a\< 1/167; \g'(0) - b\ < e = 1/108 

The proof of (i) is a simple calculation based on the estimates in Proposition [5] The integrals 
are estimated crudely, by placing an absolute value on all terms and multiplying with the length 
of the interval, 17/10. 

For (ii) we note that 

ex < pfi|| (2\\A\\\\8\\^ + \\B 1 \\S + 6\\SW 2 +\\R\\) 

(iii) At t = (s = 1.7) we have 

|ff(0) + a\ < |so(0) + o| + max KJj (0) + a 2 J 2 (0)| + £1 < 1/167 

\ a i\< a i >J=1)2 

where we used (ii), (|123j) . (|125j) and ()131j) : g'(0) is estimated in a very similar way. 
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6. The Maclaurin series of y t 

Proposition 19. The Maclaurin series of the function t i— > g(t) converges in a disk of radius at 
least Rq = 1.85. 

Proof. We define c\, i = 0, 1, 2, ... to be the Taylor coefficients of g at t = and note that 

(136) c 2 = 3c 2 , > 

The recurrence relation for the Taylor coefficients of the solution of (|121|) with initial condition 
g(0) = c = -a, g'(0) = ci = b is 

fc 

(137) (fc + l)(fc + 2)c fc+2 = 6^ Cj -c fe _ j; k> 1 

i=o 

We now take the full range of initial conditions compatible with the error range (|135|) (for simplicity 
we take the larger of the two bounds, e): cq = —a + ea\, 6\ = b + eoi where <7i, o~2 £ I = [ — 1,1], 
and calculate the formal series solution at zero for these initial data; we denote by Cj the Taylor 
coefficients thus calculated; we have 

(138) c = -a + eai, c\ = b + £0-2,0% = 3(—a + eai) 2 , £3 = 2(— a + eai)(b + ea2) + \ 

The coefficients — cq, c±, C2, 63 in fjl38|) can clearly be maximized/minimized by elementary means 
as functions of (01, 02) € I 2 . The result is 

(139) < -co < 1/5, < ci < 6/19, < c 2 < 1/8, < c 3 < 1/15 
We write 

k 

(140) (fc + 2)(fc + l)|c fc+2 | <6^|c j ||c fc _ j |; k^2 

1=0 

where Cj, i = 0, 3 are taken to be the upper bounds in (|139[) . We check that for k = 0, 1, 2, 3 we 
have 

37 

(141) |cfe| < (fc + l)/<+ 2 ; where i? = — = 1.85 
Lemma 20. The inequality (|141[) is satisfied for all fc ^ 0. 

Proof. We note that for any p > the sequence afe = (fc+ l)p k+2 , k = 0,1, 2, ... is an exact solution 
of the recurrence 

k 

(142) (fc + l)(fc + 2)a fc+2 = 6^2aja k -j,y k ^ 

3=0 

The rest is straightforward induction. 
Proposition Q1J] now follows from (|141j) . 



7. End of proof of Theorem [T] 

We have already shown that in each domain f2/, f2iUf22 the unique solution we obtained equals 
the tritronquee ht{x). After changes of variables, by matching at z — zq — 1.7e OT / 5 (i.e., at t = 
to = —1.7), we then proved that yt(z) has a convergent Maclaurin series with radius of convergence 
^ 1.85. Therefore, it follows that yt(z) is pole-free in the domain { — ^ argz ^ J, \z\ ^ 1.7} U 
{z : \z\ ^ 1.85}. By the symmetry of the solution, see Remark [3] and the Schwartz reflection 
principle, regularity also follows for argz <G [n/5,n}. 

In the process, we determined yt{0) and y' t {0) to within < 10~ 2 rigorous error bounds (see (|135[) 
and (|12ip 'l in agreement with [5]. 
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8. Appendix: The concrete expressions of various terms 



— 161 S 4 35 g 6 . _ 995 S 3 301 g 5 11 g 7 
{ ' T5 ^> ~ 64 C 2 + 1728 C 4 41472 C 6 ' r ~~ 2304 C 3 + 20736 C 5 124416 C 



7 



392 5551 S 2 1417 S 4 289 g 6 23 g 8 



(144) r 7 (C) = h 



625 9216 C 2 165888 C 4 248832 C 6 2985984 £ 



225g 2051 S 3 241 g 5 23 g 7 5 g 9 

*■ ^ rs '^ - 512C ~ 9216 C 5 " ~ 55296 C 3 " + 186624 C 7 ~ 8957952 <~ 



81 S 2 43 g 4 947 g 6 215 5 s 25 g 10 

^ ^ r9 ^ ~~ 32768 C 2 + 16384 C 4 ~ 1327104 C 6 + 23887872 C 8 ~ 859963392 C 10 

( ") 95(C) - 384 £2 + g64 ^4 " 6912 ^6 ' *(C) - 1152 £3 + 10 368 C 5 _ 124416 C 7 

95 g 3817 S 3 277 S 5 23 g 7 
*■ ^ 57 ^ ~~ 96 C 2 165888 C 4 + 41472 C 6 373248 C 8 

, A , 621g 2 1591 g 4 623 g 6 5 g 8 



1024C 3 82944 C 5 746496 C 7 995328 C 9 
15 g 3 3515 g 5 1675 g 7 125 g 9 



^ 15 °' ) 99 ^ 2048 C 4 884736 C 6 + 23887872 C 8 429981696 C 10 

, . __269 5 3 61 g 5 1691 g 4 353 g 6 

( j 5(Cj ~ 576 C + 10368 C 3 ' { ' ~ 20736 C 2 + 497664 C 4 

95 g 3 1915 g 5 25 g 7 25 g 4 125 g 6 625 g 

( ) -^HO C-7C /" _ O/IOOOQ AM 0700/10 AS ' -^ 8 (C) r. co ,9 ~~ Ifirooo AA • 



576 C 248832 C 3 373248 C 5 ' 768 C 2 165888 C 4 143327232 C 

(153) Eu{p) = f «!! + « ^ p - 3 / 2 + f «E» + « 
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288 15552 ^ 1 \ 20736 995328 

/ 95|g| 3 1915|g| 5 5|g| 7 \ _ 6/2 / 25[g| 4 125|g| 6 625]5| L \ _ 3 
V 288 373248 186624 / P V 768 331776 429981696 1 



(154) M q = Y, M ^l- 



1U 9J/ 

where 

(155) M 9 , 3 = J| |g| + || ,g, 3 + ^ ISf , M qA = i±W 2 + ^^ + dk^ 
^ M - = I 1*1 + 82H0 ^ + 20^0 I 5 ' 5 + 18^40 I 5 ' 7 

(158) M,, 7 = |g| 3 + |g| J + 14 33 6 2 7 7 5 232 is 1 ! 7 + ag^ggo^g l^l 9 

7 

(159) M L ,, = ^M L , g , jp - j/2 , 
where 

( 16 °) M ^ = M + SES + sk ^ < M ^ = SI ^ + 41I2 ^ + 497^64 '^ 

< 161 > M - - § I 5 ' + tUs I 5 ' 3 + iSi ^ 5 + Telle 
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^ M ^ = S 15,2 + So ^ + srilo ^ + 995W ' 5 ' 8 

15 3515 1675 125 

(163) Ml, 9 , 7 = \S\ 3 + 486604g |5| 5 + 131383296 l^l 7 + 2364899328 ^ 

2 

(164) M G ,i = ^ mj -,ip- J ' /2 

(165) mo, = ™ ( 1 + V2U {S f + » |g| 8 + \Sf + . 1417 



3125 \ J 11520 1 1 3732480 1 1 311040 1 1 207360 1 1 

(166) mil = 2L\s\ + |sf + 2*L \ S \ 5 + — 5 - — \s\ 9 + \s\ 7 

y ' ' 256 1 1 13824 1 1 82944 1 1 13436928 1 1 279936 1 1 

947 . |6 25 ,io 215 . |8 43 , |4 81 , |2 

(167) roa ' 1 - 2322432 |S| + 1504935936 |S] + 41803776 ]S] + 28672 |5] + 57344 |S| 

(168) Mg2= 16^ + ^ + ^ + /9^ 11^ 30^\ 1/2 
v ; ' 4320 160 20736 \ 6912 373248 62208 J P 

(1 69) Mo ,, _ ^ |sr + ^| S |« + 52| SP+ (^ |S |. + ^ |S| . + _J£_| Sr )^ 

161 g 4 53 g 2 35 g 6 . , 23 g 7 35 g 5 1631 g 3 

^ ^ 5 ^" 1728 C 5 64 C 3 41472 C 7 ' 62208 C 8 + 768 C 6 2304 C 4 

11 g 8 301 g 6 995 g 4 
*■ ^ 7 ^ ~ 373248 C 9 + 62208 C 7 6912 C 5 

161 g^ _ _53_ _ 35 _ 47 g^ _ 1631 _913_ g^ 

(1"2) U 6 (Q - 345g ^ 3 12g ^ g2944 ^ 5 , "6(C) - 12441 g ^ 6 4U72 £ 4 + 460g ^ 

173 g 8 3479 g 6 1843 g 4 11 g 9 301 g 7 995 g 5 

( 73) U 7 (g - 7464g6 ^ 7 ^ 5 + 460g ^ , "8(C) - 559g72 ^ 8 g3312 ^ 6 + 1Q368 ^4 

161g 4 53g 2 5g 6 _ 23g 7 35g 5 1631g 3 

' T5 ^' ~ 8640C 5 + 192C 3 + 41472C 7 ' ^ ' ~ 497664C 8 ~ 4608C 6 + 9216C 4 

,„ rl llg 8 43g 6 199g 4 

(!75) MO = oo.nooo.O " m , n . W + 



3359232C 9 62208C 7 6912C 5 



.3 



- _ 161 g 4 265 gj^ 25 g 6 - _ 23 g 7 35 g 5 1631 S_ 

( ' 7[ ^' ~ 3456 C 5 + 384 C 3 + 82944 C 7 ' 165888 C 8 1536 C 6 + 3072 C 4 

- 77 g 8 301 g 6 1393 g 4 
( ' 9 ^> ~~ 6718464 C 9 ~ 124416 C 7 + 13824 C 5 

161 g 4 _53_g2_ 7 g 6 47 g 7 1631 g 5 913 g 3 

1/6 ^ " " 10368 C 3 + 128 C + 82944 C 5 ' ~ ~ 746496 C 6 + 165888 C 4 ~ 9216 C 2 

173 g 8 3479 g 6 1843 g 4 11 g 9 301 g 7 995 g 5 

( ) - 5225472 £7 + 622080 C 5 13824 C 3 ' ~ 4478976 C 8 + 559872 C 6 41472 C 4 

~ tn- 805 -S" 1 265 35 g 6 - rn _ 47 g 7 1631 g 5 913 g 3 

(180) w(C) - "20736 C 3 + 256 C + 165888 C 5 ' ~ "248832 C 6 + 55296 C 4 ~ 3072 C 2 

(IKU ~ 173 g 8 24353 g 6 12901 g 4 . 11 g 9 301 g 7 995 g 5 

( ) "9(C) - - 1492992 ^7 + 1244160 C 5 ~ 27648 C 3 ' ~~ _ 1119744 C 8 + 139968 C 6 ~ 10368 C 4 

161 g^ _53_ Ef_ 1 g^ _ 1001 g^ _ 65 g^ _ 17 g^ 

*■ ' Pb ^' ~ 25920 C 3 + 192 C + 41472 C 5 '' ~ 829440 C 4 18432 C 2 2985984 C 6 

_ 17 g 6 185 g 8 137 g 4 199 g 5 43 g 7 11 g 9 

( ) P7(C) - lg440 ^ 5 - 470292 48 £7 ~ 4608 C 3 ' " "41472 C 4 + 559872 C 6 ~ 40310784 C 8 
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(184) 

== 15,4 + M ^ + Ilk & + (^84 & + siSo ' 5 ' 5 + diS |S|3 ) ^ 



V 47029248 1 1 4608 1 1 19440 1 ' }' \ 559872 1 1 40310784 1 1 41472 



P" 1/2 



(185) M - 15,6 I 161 ' S I 4 + + /" 2401 l g ! 5 + 761 l^| 3 , 4 97|g| 7 

1 ' 6912 6480 96 V 92160 2048 2985984 

/miJSf 3083|S| 8 3673|g| 4 \ _i /199jSf 559|S| 7 187|g| 9 \ _ 3/2 

+ V 155520 + 47029248 + 13824 J P + \ 4608 + 559872 + 40310784 J 9 

< I86 > M " - shs |S| * + 55^68 > s i* + + tH |S|2 + 1 

+ (l^o |s| ° + liS5 |if + H |S| + itkwn |S| ° + mIm |s| ') 

, / 2 5 c,,io , 2 15 |^,8 947 | , 6 43 , , 4 27 , „, 2 \ -i 

V 2579890176 1 1 71663616 1 1 3981312 1 1 49152 1 1 32768 11 J P 

m Jfc ,« . Jf, isr + W + i§ | 4 f + ( -£| W + ^ W + w) 
<"»> M - " mk |S| ' + SS5S' S " , + * S <" + (S |S| * + ISH |S|S + rJSso i s i 7 ) 

/ 12139 4 2623 6 671 8 \ , 

\ 290304 1 1 2612736 1 1 141087744 1 1 ) P 

8.1. Values of intermediate constants for p = 3. The numerical values of these constants might be 
helpful to the reader who would like to double-check the estimates. These are: 

J M = 0.282580 ,j m = 0.64374..., Y 1>M = 1.16314..., Y hRiM = 0.132618- , E M = 0.0490292-. 
Z2,r,m = 0.54226 , z 2 , A/ = 0.91863 ,M 9 = 0.066702 .., M L , q = 0.075708 ■ -, V M = 0.2239 ,T M = 0.0385- 
Mi = 1.13838 ,M 2 = 0.04303 -,M 3 = 0.28346 ,M 4 = 0.45227 M 5 = 0.05430 , M 6 = 0.00231 -,M 7 = 0.02018 
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